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We calculate the Cherenkov radiation from an e + e~ pair 
at small separations, as occurs shortly after a pair conversion. 
The radiation is reduced (compared to that from two indepen- 
dent particles) when the pair separation is smaller than the 
wavelength of the emitted light. We estimate the reduction in 
light in large electromagnetic showers, and discuss the impli- 
cations for detectors that observe Cherenkov radiation from 
showers in the Earth's atmosphere, as well as in oceans and 
Antarctic ice. 



I. INTRODUCTION 

Cherenkov radiation from relativistic particles has 
been known for over 70 years [1]. However, to date, al- 
most all studies have concentrated on the radiation from 
individual particles. Frank [2], Eidman [3] and Balazs 
[4] considered the Cherenkov radiation from electric and 
magnetic dipoles, but only in the limit of vanishing sep- 
arations d. Their work was nicely reviewed by Jelley [5] . 

Several more recent calculations have considered 
Cherenkov radiation from entire electromagnetic show- 
ers, in the coherent or almost coherent limit [6]. The 
fields from the e + and e~ largely cancel, and the bulk of 
the coherent radiation is due to the net excess of over 
e + (the Askaryan effect) [7]. Hadronic showers produce 
radiation through the same mechanism [8] . Coherent ra- 
diation occurs when the wavelength of the radiation is 
large compared to the radial extent of the shower; for 
real materials, this only occurs for radio waves. 

Here, we consider another case, the reduction of ra- 
diation from slightly-separated oppositely-charged co- 
moving pairs. This includes e + e~ pairs produced by 
photon conversion. When high-energy photons convert 
to e + e~ pairs, the pair opening angle is small and the e + 
and e~ separate slowly. 

Near the pair, the electric and magnetic fields from the 
e + and e~ must be considered separately. However, for 
an observer far away from the pair (compared to the pair 
separation d), the electric and magnetic fields from the e + 
and e~ largely cancel. Cherenkov radiation is produced 
at a distance of the order of the photon wavelength A 
from the charged particle trajectory. So, for d < A, can- 
cellation reduces the Cherenkov radiation from a pair to 
below that for two independent particles. For a typical 
pair opening angle m/k, where k is the photon energy 
and m the electron mass, without multiple scattering, 
A > d for a distance kA/m. For blue light (A = 400 
nm) from a 1 TeV pair, the radiation is reduced until 



the pair travels a distance of 40 cm (neglecting multiple 
scattering). 

A similar cancellation effect was observed for energetic 
(~ 100 GeV) e + e~ pairs in nuclear emulsions [9]. Ion- 
ization from newly created e + e~ pairs is reduced when 
the pair separation is less than the screening distance for 
ionization in the target. 

In this paper, we calculate the Cherenkov radiation 
from e + e~ pairs, simulate optical radiation from pairs 
follow realistic trajectories, and consider the radiation 
from electromagnetic showers. We consider two classes 
of experiments: underwater/in-ice neutrino observatories 
and air Cherenkov telescopes. 



II. CHERENKOV RADIATION FROM PAIRS 

Cherenkov radiation from closely spaced e + e~ pairs 
can be derived by extending the derivation for point 
charges, by replacing a point charge with an oppositely 
charged, separated pair. We sketch the derivation for 
radiation from point charges, review previous work on 
radiation from infinitesimal dipoles, and derive the ex- 
pression for Cherenkov radiation from a closely-spaced 
co-moving pair. 

We follow the notation and derivation from Ref. [10]. 
In Fourier space, the charge density p and current density 
J from a point charge ze propagating with speed v in the 
x\ direction can be written as 



P(k,uj) = —S(lj - kiv) 
J(k, u>) = vp(k, uj) 



(1) 



where k is the wave vector and ui the photon energy. This 
current deposits energy into the medium through elec- 
tromagnetic interactions. We use Maxwell's equations 
beyond a radius a around the particle track, where a is 
comparable to the average atomic separation. Then, by 
conservation of energy, the Cherenkov radiation power is 
equal to the the energy flow through a cylinder of this 
radius, giving 



dE 
dx 



-caRe 



B^E^dw 



(2) 



Ei is the component of E parallel to the particle track, 
and B3 is the component of B in the X3 direction, eval- 
uated at an impact parameter b at a point with X2 = b, 
X3 = 0. We omit the time-phase factors for brevity. 
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Using the wave equations in a dielectric medium and 
the definition of fields, then integrating over momenta, 
which eliminates the space-phase factors, one finds 



where 
Similarly, 



K (Xb) (3) 



ze (2 

v \tt J e(ui) 



B 3 (u>) = e(w)0E 2 (u) . 



(4) 



K and K\ are the zeroth and first order modified Bessel 
functions of the second kind. 

The far-field radiation depends on the asymptotic form 
of the energy deposition at |Ao| ^> 1. For > l/^/e(w) 
for real e(o>), A is completely imaginary. The asymptotic 
contribution of the Bessel functions in the integrand of 
dE/dx is finite, giving the well-known expression for the 
Chcrenkov radiation 



dE 
dx 



{zef 



J 

Je(u) 



X / LU 
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(5) 



duJ . 



Note how a has dropped out [10, Ch. 13]. The derivation 
of this Cherenkov radiation may be expanded to give the 
field from a pair. 

The radiation from an e+e~ pairs depends on two pa- 
rameters: the separation d and the angle between the 
direction of motion and the orientation of the pair. For 
relativistic pairs created by photon conversion, the trans- 
verse (to the direction of motion) separation is important; 
the longitudinal separation of a highly relativistic pair 
can be neglected, due to Lorcntz length contraction. 

Balazs [4] provided an expression for Cherenkov radia- 
tion from an infinitesimal dipole D oriented transverse to 
its momentum. The fields are approximated by by a lin- 
ear Taylor expansion of the corresponding point-charge 
fields: 



where d is the effective pair separation, so D = zed. 
Then, following the same steps as in the point-charge 
case, Balazs finds 



dE 
dx 
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lc( U )>i/tp \ P £ M 



(6) 

dto . 



For a point dipole oriented parallel to its direction of 
motion, the radiation is negligible for < 1 [5]. 

To compute the Cherenkov radiation for finite separa- 
tions d, let us consider a pair moving in the +x direction. 
The pair lies entirely in the transverse plane y-z, with the 
line between them making an angle a with respect to the 
y-axis. Then, generalizing Eq. (1), the charge density 
from the pair is 



p(k,u) = ^6{u-k 1 v)[e-^ k2y +- k3Z +^-e- i ^ y — k3Z ^] 



The two charges have positions relative to the center of 
mass 



d 



d 



V+=2 COSa z + = --sma 



y_ = — cos a z_ — — sin a . 
y 2 2 

The angle a is the relative azimuth between the line con- 
necting the two charges and the azimuth of observation. 
The generalization of E\(u) of Eq. 3 is 



£i(w) 



-0' 



—izeuj 
x \K (\b-) - K (\b+)} 



1 



(7) 



where 



b± = y — sin a + (b ± — cos a) z . 

As before, we take |Aa| 3> 1 and a < b, so we need only 
consider d <C b; there is little interference for d > b. 
Therefore, we can simplify using 

b± ~ b ± ^ cos a . 

Then, as before, considering completely imaginary A and 
|Ao| > 1, 




^|A| cos a 



(8) 



and a similar expression for ^3(0;). Here we have taken 
6_ ~ b + ~ b in the denominator. 

At a = ±7r/2, E\(w) = 0. The Cherenkov radiation is 
no longer symmetric about the direction of motion, and 
vanishes at right angles to the direction of the dipole. As 
the charge separation increases (or the wavelength de- 
creases), the angular distribution evolves from two wide 
lobes into a many-lobed structure, as shown in Fig. 1. 
After integration over even a narrow range of ui or d, the 
angular distribution becomes an almost-complete disk, 
with two narrow zeroes remaining at a direction perpen- 
dicular to the dipole vector. 
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FIG. 1. The azimuthal angular distribution (transverse to 
the direction of motion) of Cherenkov radiation for 500 nm 
photons from a pair of charges oriented as shown in the Fig- 
ure. Distributions are shown for pair separations 100 nm 
(solid line), 1 fim (dashed line) and 5 /im (dotted line), with 
sj — n — 1.3 and ft = 1. 

After assembling the pieces, and averaging over a, we 
find the generalization of Eq. (5), 



dE\ {zef 



dx J 



uj 1 



M>l//3 

1 



x 2 [1 - J (Ad)] 



(9) 



Here Jo is the zeroth order Bessel function of the first 
kind. For Xd <C 1, this reproduces Eq. (6). For Xd^> 1, 
the dE/dx is twice that expected for an independent par- 
ticle (Eq. (5)). The transition is shown in Fig. 2. As the 
emission wavelength A approaches d, the pair spectrum 
converges to the point-charge spectrum in an oscillatory 
fashion, characteristic of the Bessel function. For certain 
values of Xd, the radiation exceeds that of two indepen- 
dent charged particles. 

For th e rem ainder of the paper, we assume that media 
satisfy y e ( UJ ) = n , where n is independent of frequency. 
In realistic detection media, any variation of n with fre- 
quency is small, and would have little effect on Cherenkov 
radiation from relativistic particles. 

With real e + e~ pairs, two effects should be considered. 
Electromagnetic radiation is not emitted instantaneously, 
but occurs while the radiating particles travel a distance 
known as the formation length, If. For Cherenkov radi- 
ation, l f = A/sm 2 (6» c ) = Ae(3 2 /{e(3 2 - 1) [11], depends 
only on the Cherenkov emission angle, 9c and the pho- 
ton wavelength; If depends only slightly (through [3) on 
the electron energy. 
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FIG. 2. The spectrum of Cherenkov radiation at /3 = 1, 

sj e(uj) = n = 1.3. Solid line is for e + e~ with the particles 

considered independently, and the dashed lines are for pairs 

treated coherently, with separations 100 nm, 1 /jm and 5 /im. 

While the pair is covering the distance If, the pair sep- 
aration will change by an amount Ad = / f sin (9) , where 
9 is the angle between the e + and e~ velocity vectors. 
Since 9 is of order 1/7, Ad/d <C 1, so the change in 
separation is not significant. 

Second, the Cherenkov radiation produced at a point 
(x-coordinate) depends on the fields emitted by the 
charged particles at earlier times, when d may be dif- 
ferent than at the point of radiation. For full rigor, these 
retarded separations should be used in the calculation. 
Again, this has a negligible effect on the results. 



III. RADIATION FROM e+e" 
SHOWERS 



PAIRS IN 



Many experiments study Cherenkov radiation from 
large electromagnetic showers. The radiation from a 
shower may be less than would be expected if every parti- 
cle were treated as independent. We use a simple simula- 
tion to consider 300 to 800 nm radiation from electromag- 
netic showers. This frequency range is typical for pho- 
tomultiplier based Cherenkov detectors; at longer wave- 
length, there is little radiation, while shorter wavelength 
light is absorbed by the glass in the phototube. 

We simulated 1000 7 conversions to e + e~ pairs with 
total energies from 10 8 to 10 20 eV. Pairs were pro- 
duced with the energy partitioned between the e + and 
e~ following the Bethe-Heitler differential cross section 
da ps E±(l — E±), where E± is the electron (or positron 
energy) [13]. At high energies in dense media (above 10 16 
eV in water or ice), the LPM effect becomes important, 
and more asymmetric pairs predominate [12]. The pairs 
are generated with initial opening angle otm/k; the fixed 
angle is a simplification, but the pair separation is domi- 
nated by multiple scattering, so it has little effect on our 
results. 

The e~ and e + are tracked through a water medium 



(with n = y/e = 1.3) in steps of 0.02Ao, where Xq is the 
radiation length, 36.1 cm in water. At each step, the par- 
ticles multiple-scatter, following a Gaussian approxima- 
tion [14, Ch. 27]. The particles radiate bremsstrahlung 
photons, using a simplified model where photon emission 
follows a Poisson distribution, with mean free path Xq. 
Although this model has almost no soft bremsstrahlung, 
soft emission has little effect on Cherenkov radiation, 
since the electron or positron velocity is only slightly af- 
fected. 

At each step, we compute the Cherenkov radiation for 
each pair. They are treated coherently when d < 2A; at 
larger separations the particles radiate independently. 

As shown in Fig. 3, the particles in lower energy pairs 
(< 10 10 eV) radiate almost independently. In contrast, 
the radiation from very high energy pairs (> 10 15 eV) is 
largely suppressed. The broad excursions slightly above 
unity occur when Jo (Ad) > 1 for many of the scattered 
pairs. 
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FIG. 3. Average Cherenkov photon emission rate for pairs 
with energies from 10 10 (leftmost, dot-dashed curve) to 10 15 
eV (rightmost, dotted curve) vs. the distance travelled by 
the pair in water, relative to emission from two independent 
particles. 



IV. IMPLICATIONS FOR EXPERIMENTS 

At least two types of astrophysical observatories de- 
pend on Cherenkov radiation. Water and ice based neu- 
trino observatories observe Cherenkov radiation from the 
charged particles produced in neutrino interactions, and 
air Cherenkov telescopes look for 7-ray induced electro- 
magnetic showers in the Earth's atmosphere. 

Current neutrino observatories can search for electron 
neutrinos with energies above 50 TeV (for v^, the thresh- 
old is much lower) [16]. They use large arrays of photo- 
multiplier tubes to observe the Cherenkov radiation from 
v e induced showers. For water, n ~ 1.3, Fig. 3 shows 
that Xd < 1 while the pair travels significant distances. 
Ice is similar to water, with a slightly lower density; n of 
ice depends on its structure, and is typically « 1.29 [18]. 



To quantify the effect of Cherenkov radiation from u e 
interactions, we use a toy model of an electromagnetic 
shower. The shower evolves through generations, with 
each generation having twice as many particles as the 
preceding generation, with half the energy. Each gen- 
eration evolves over a distance of Xq] other simulations 
have evolved generations over a shorter distance (In 2)Xo, 
leading to a more compact shower [17]. In these showers, 
most of the particles are produced in the last radiation 
lengths. 

Fig. 4 shows the Cherenkov radiation expected from a 
model 10 20 eV shower with coherent Cherenkov radiation 
(solid line) and in a model where all particles radiate in- 
dependently (dotted line). This model does not include 
the LPM effect, so it should be considered only illustra- 
tive. The LPM effect lengthens the high-energy (above 
a few 10 15 eV) portion of the shower. By spreading the 
shower longitudinally, the LPM effect will give the elec- 
trons and positrons more time to separate, and so will 
somewhat lessen the difference between the two results. 
However, it is clear from Fig. 4 that coherence has a 
significant effect for the first w 22 generations. Since the 
front of the shower contains relatively few particles, it will 
not affect the measured energy; the change in number of 
radiated photons (and hence on the energy measurement) 
should be less than 1%. 




35 40 
Depth [XO] 

FIG. 4. Cherenkov radiation from a 10 20 eV shower in 
water, using the Heitler toy model, versus shower depth 
(smoothed). The two curves compare the radiation for e + e~ 
calculated as independent particles and as coherent pairs. 

However, the suppression will affect the apparent 
length of the shower. For the first w 8 generations, the 
shower will emit less light than a single charged parti- 
cle. Because of the LPM effect, each of these genera- 
tions (with mean particle energy E g above a few 10 15 
eV) develop over a distance X = X n yjE g /bEi,PM, where 
Elpm = 278 TeV for water is the effective LPM energy 
[17], greatly elongating the shower. So, the first 8 gen- 
erations include most of the length of the shower. So, 
the suppression of Cherenkov radiation hides the initial 
shower development, making the shower appear consid- 
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erably more compact. The reduction in early-stage radi- 
ation should help in separating electron cascades from 
muon-related backgrounds, especially muons that un- 
dergo hard interactions, and lose a large fraction of their 
energy. 

Atmospheric Cherenkov telescopes like the Whipple 
observatory study astrophysical 7-rays with energies 
from 100 GeV to 10 TeV. These telescopes observe 
Cherenkov radiation from pairs in the upper atmosphere; 
for a 1 TeV shower, the maximum particle density occurs 
at an altitude of 8 km above sea level (asl) [15], where the 
density is about 1/3 that at sea level. Since n— 1 depends 
linearly on the density, at 8 km asl n—1 f=s 1 x 10~ 4 , so for 
500 nm photons radiated from ultra-relativistic particles, 
Ac? < 1 only for d < 6 /um. In this low-density medium, 
the effect of the pair opening angle is significant and mul- 
tiple scattering is less important. Pairs with k < 1 TeV 
will separate by 30 /<m in a distance less than 30 meters; 
at 8 km asl, this is 3% of a radiation length. This dis- 
tance is too short to affect the radiation pattern from the 
shower. 

Cherenkov radiation is also used in lead-glass block 
calorimetry, and in Cherenkov counters for particle iden- 
tification; their response to photon conversions may be 
affected by this coherence. 

Although the reactions are slightly different, a similar 
analysis applies to the reduction of ionization by e+e~ 
pairs. Perkins observed that the ionization from pairs 
with mean energy 180 GeV in emulsion was surppressed 
for the first w 250 /um after the pairs were created [9]. 
With Xq = 3 cm (typical for emulsion), the e + and e~ 
trajectories will be about 4 nm apart after travelling 250 
/an. For relativistic particles, the screening distance (ef- 
fective range for dE/dx) is determined by the plasma 
frequency of the medium, lo p . For silver bromide, the 
dominant component of emulsion, fno v = 48 cV [10] (in 
a complete emulsion, fraj p will be slightly lower). This 
yields a screening distance c/uj p = 4 nm, which is very 
close to the calculated separation. 
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V. CONCLUSION 



We have calculated the Cherenkov radiation from e + e~ 
pairs as a function of the pair separation d. When d 2 < 
v 2 1 (uj 2 [1 — P 2 e(co)]), the radiation is suppressed compared 
to that from two independent particles. 

This suppression affects the radiation from electromag- 
netic showers in dense media. Although the total radi- 
ation from a shower is not affected, emission from the 
front part of the shower is greatly reduced; this will af- 
fect studies of the shower development, and may affect 
measurements of the position of the shower. 
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